Theoretical investigation of polarization-compensated II-IV/I-V perovskite 

superlattices 
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Recent work suggested that head-to-head and tail-to-tail domain walls could be induced to form 
in ferroelectric superlattices by introducing compensating "delta doping" layers via chemical substi- 
tution in specified atomic planes [Phys. Rev. B 73, 020103(R), 2006]. Here we investigate a variation 
of this approach in which superlattices are formed of alternately stacked groups of II-IV and I-V 
perovskite layers, and the "polar discontinuity" at the II-IV/I-V interface effectively provides the 
delta-doping layer. Using first-principles calculations on SrTiOs/KNbOs as a model system, we 
show that this strategy allows for the growth of a superlattice with stable polarized regions and 
large polarization discontinuities at the internal interfaces. We also generalize a Wannier-based def- 
inition of layer polarizations in perovskite superlattices [Phys. Rev. Lett. 97, 107602 (2006)] to the 
case in which some (e.g., KO or Nb02) layers are non-neutral, and apply this method to quantify 
the local variations in polarization in the proposed SrTiOa/KNbOa superlattice system. 

PACS numbers: 77.80.Dj 77.22.Ej 77.84.Dy 81.05.Zx 
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Ferroelectric materials have been the subject of in- 
creasing theoretical and experimental study in recent 
years. In particular, multicomponent superlattices 
based on ABO3 perovskites have been shown to pos- 
sess many interesting properties (see Ref. [l| and refer- 
ences therein). Currently, heterointerfaces between dif- 
ferent ABO3 perovskites are also the subject of intense 
investigationi^iii^i^ Given the huge number of possible 
superlattice arrangements, it is prudent to investigate 
systems of potential interest using first-principles calcu- 
lations to confirm a given system possesses the desired 
properties before experimental work is performed.'^ 

In recent work, Wu and Vanderbilt^ introduced a novel 
concept in which 180° head-to-head (HH) and tail-to-tail 
(TT) domain walls are induced to form in a ferroelec- 
tric superlattice via the insertion of compensating "delta 
doping" layers. For example, in a II-IV ABO3 perovskite, 
column III or V ions could replace the IV ions in one BO2 
layer, inducing the formation of a HH or TT domain 
wall respectively. Usually the formation of a HH or TT 
domain wall perpendicular to the polarization direction 
would entail an unacceptable Coulomb energy cost, or 
cause the domain wall to become metallic in order that 
free carriers could compensate the domain wall. How- 
ever, it was shown that the delta-doping layers could be 
arranged to compensate the polarization bound charge 
and allow a structure in which the ferroelectric domains 
are polarized in opposite directions and are separated by 
HH and TT domain walls. 

Here we examine the possibility of forming a multicom- 
ponent perovskite superlattice with similarly large dis- 
continuities in the local electric polarization by making 
use of alternating II-IV and I-V perovskite constituents. 
In this case, the "polar discontinuity" associated with 
the II-IV/TV interface plays the role of delta-doping 
layer and compensates the polarization bound charge at 
the interface. (Clearly a similar strategy could be ap- 
plied to II-IV/IITIII perovskite superlattices.) The re- 



sulting structure is not switchable, but is locally polar- 
ized and has strongly broken inversion symmetry. We 
demonstrate this concept via first-principles calculations 
on SrTiOa/KNbOs as a prototypical system, showing 
successful compensation and robust formation of locally 
polarized regions. Superlattices of this type are shown 
to remain insulating to rather large layer thicknesses. 
Finally, we clarify how the Wannier-based definition of 
layer polarization in perovskite superlattices introduced 
in Ref. [l^ can be generalized for a system having non- 
neutral AO or BO2 constituent layers, and apply this to 
the SrTiOa/KNbOa system to map out the local varia- 
tions in polarization. 

The type of superlattice structure we have in mind is 
illustrated in Fig. [T] Here, two unit cells of KNbOs (in 
the sequence Nb02-KO-Nb02-KO) repeatedly alternate 
with two unit cells of SrTiOa (in the sequence Ti02-SrO- 
Ti02-SrO) during growth. While each added KNbOa or 
SrTiOa unit cell is neutral, in KNbOs this neutrality re- 
sults from the cancellation of charges on the K0~ and 
NbOj layers, while in SrTiOs the individual layers are 
neutral. When the layers are assembled as in Fig. [2 
the presence of the "polar discontinuity" introduces ef- 
fective compositional charges of ±1/2 at the Nb02/SrO 
and Ti02/K0 interfaces respectively, as shown. Simi- 
lar effects have recently been extensively discussed for 
SrTiOa/LaAlOa and related interfaces i^i^^i^i^ Intuitively, 
we may regard each K0~ layer in bulk KNbOs as being 
half-compensated from each of its two immediate NbOj 
neighbors, whereas at the Ti02/K0 interface the KO^ is 
only half-compensated because it has only a single NbO J 
neighbor. The resulting compositional charge densities 
are (Tcomp = ±e/2a'^ at the Nb02/SrO and Ti02/K0 
interfaces respectively. 

Two possibilities for cancelling these compositional 
charges in a I-V/ITIV superlattice are shown in Fig. [2l 
In Fig.[2lja) we assume that both I-V and II-IV materials 
are ferroelectrics, and to take an extreme case we assume 
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FIG. 1: (Color online) Sketch of a possible II-IV (SrTiOs) / 
I-V (KNbOs) superlattice grown along [001]. 



their spontaneous polarizations are equal {Ps^^ = Ps^'') 
at the specified in-plane lattice constant. Arranging the 
polarizations so that they alternate up and down as 
shown, the ideal compensation (Tbound + Ccomp = be- 
tween polarization bound charge and compositional in- 
terface charge is realized when P^^ = Ps^'^ = e/4a^. A 
second scenario, shown in Fig. [2Ib), would result from 
the alternation of I-V and II-IV materials, only one of 
which is ferroelectric, while the other is paraelectric. In 
this case, ideal compensation requires Pg = ejla^ for the 
ferroelectric component. Intermediate cases, with the al- 
ternation of a strong and a weak ferroelectric, are also 
possible. 

In the remainder of this paper, we focus on the scenario 
sketched in Fig. H^b) as reahzed in the SrTiOa/KNbOa 
superlattice system, using density-functional calculations 
to demonstrate the compensation mechanism proposed 
above. The choice of SrTiOa as the II-IV paraelectric 
component is motivated by the fact that it is a well- 
studied materiaUiii^ and by the common use of SrTiOs 
as a substrate for growth of thin perovskite films. With 
this in mind, we assume coherent epitaxy of our super- 
lattices on SrTiOa, so that the in-plane lattice constants 
are constrained to the experimental value ao = 3.905 A 
of bulk SrTiOa. As explained above, ideal compensa- 
tion would require that the ferroelectric I-V component 



should have a spontaneous polarization of Pg = e/2aQ — 
0.525 Cjvc? . We have chosen KNbOa for the I-V compo- 
nent because it provides a reasonable match to this value. 
While bulk KNbOa is a rhombohedral ferroelectric with 
polarization along (111) at low temperature, its lattice 
constant is nearly 3% larger than that of SrTiOs. The 
calculations of Dieguez et al}^ have shown that SrTiOs 
should become a tetragonal ferroelectric with polariza- 
tion along (001) when compressed in-plane to fit to the 
SrTiOs lattice constant. Moreover, those same calcula- 
tions indicated that the polarization of KNbOa would 
be ~0.45C/m^ under these conditions, which is fairly 
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FIG. 2: (Color online) Two possible II-IV/I-V superlattice 
arrangements yielding compensating heterointerfaces and sta- 
bilized ferroelectric discontinuities, (a) Both materials are 
ferroelectric, with antiparallel polarizations along the growth 
direction, (b) One material is paraelectric and the other 
ferroelectric. Green and red interface charges denote polar- 
discontinuity and polarization-induced charges, respectively. 



close to the target value. KNbOs also has the advan- 
tage of being a commonly used and very well-studied I-V 
perovskite J^ii^ 

The open-source plane-wave density-functional code 
pwscpii was used for the calculations, with the local- 
density approximation to exchange and correlatioi*i^ and 
use of ultrasoft pseudopotentials.^^ Because we have no 
reason to expect the appearance of in-plane polarization 
components in this system (see above) , we have assumed 
tetragonal P4mm symmetry throughout. 

We first calculate the theoretically optimized value of 
the lattice constant ao for bulk SrTiOa using an 8 x 8 x 8 
k-point grid and a plane-wave energy cutoff of 30 Ryd- 
berg, obtaining a value of value of 3.849 A. This leads to 
a theoretical ideal matching condition of Ps=0.54C/m^. 
Using the same k-point grid and cutoff to study bulk 
KNbOs in P4mm symmetry with its in-plane lattice con- 
stant constrained to this oq, we calculate its spontaneous 
polarization to be 0.42 C/m^. These results are in good 
agreement with previous workii^ 

Of primary importance is whether the compensation 
of the bound charge is sufficient to maintain the insu- 
lating nature of the system as the supercell size is in- 
creased. We consider superlattices consisting of n unit 
cells of SrTiOs alternating with n unit cells of KNbOs, 
so that the supercell contains lOn atoms. Figure [T] illus- 
trates the case of n=2. Relaxations of the multilayered 
supercells are performed for values of n ranging from 1 
to 5. The plane- wave energy cutoff is 30 Rydberg in all 
cases, and the k-point grid is 8 x 8 x M with Af = 4 for 
n = 1 and M = 2 for n > 2. 

We find that the system remains insulating in all these 
cases. The density of states for n = 2 is shown in Fig. [31 
showing a clear gap between valence and conduction band 
states. There is, however, a gradual closing of the band 
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FIG. 3: Calculated electronic density of states for a super- 
cell with n = 2 (two SrTiOa layers and two KNbOa lay- 
ers). Inset shows the the gradual reduction of band gap for 
nSrTiOa/nKNbOa supercells as n increases. 



gap as n is increased, as may be expected from the im- 
perfect charge compensation. The calculated band gap 
as a function of the supercell size is shown in the inset 
to Fig. [31 The results for n > 2 suggest a nearly linear 
decrease in band gap with increasing n. This reduction 
is rather modest; a simple linear extrapolation suggests 
that the system would not become metallic until n ~ 32. 

In order to understand the resulting local polariza- 
tions in these multilayered systems, it is useful to start 
from a simple model and test whether its predictions are 
borne out by a more detailed analysis. Due to the im- 
perfect compensation in the system, if one starts by as- 
sembling regions of spontaneously polarized KNbOa al- 
ternating with regions of unpolarized SrTiOa, one finds 
planar charge densities icr at the interfaces, with a — 
Ps — ejlc? = — 0.12C/m^. Assuming that the screen- 
ing of this charge can be treated using a linear dielectric 
analysis and that the thicknesses of the two constituents 
are approximately the same, one finds a screened charge 
density of (Tgcr — 2cr/(ei + £2), where £1^2 are the dielec- 
tric constants of KNbOs and SrTiOs respectively. The 
screened polarizations then become -PKNbOs = ^s^cr(ei — 
+ £2) and PsrTiOa = ^^{^2 - + £2)- Since 

(T < 0, we expect PKNbOa to be enhanced slightly beyond 
its spontaneous value, and PsrTiOs to have a small value 
of the opposite sign. 

To investigate the correctness of this picture, we per- 
form more precise calculations of the local polarization 
profile. An accurate method for obtaining layer polariza- 
tions was introduced in Ref. 10, where a one-dimensional 
Wannier analysis^i was employed. To do this, the over- 
lap matrices Mmn = {umk\ Un,k+b) between neighboring 
k-points along strings in the z direction are computed, 
where Umk is the periodic part of the Bloch function 
ipmk- The singular value decomposition M = VT,W^ is 
used to obtain M = UW'' , which is exactly unitary. The 
eigenvalues Am of the product A = M'^'^^ of these ma- 
trices along the k-point string yield the Wannier centers 



as Zm = (— c/27r)Im InAm- These Wannier centers form 
into "sheets" of charge that are localized in the growth 
direction but delocalized in the plane. For each layer j, 
we define the layer center zqj to be the average position 
of the ions associated with that layer. The "intralayer 
polarization" is then given by 

where S* = is the basal cell area, Qt is the core charge 
of ion T belonging to layer j, Rtz is the position of the 
ion measured relative to zqj , and Zm is the position of 
the Wannier center Zm relative to zqj after a and ky 
average over wavevector strings. 

As long as each layer is neutral, the total polarization 
of the supercell can be obtained just as (J^jPj)/'^ where 
c is the supercell lattice constant, because the dipole mo- 
ment of a neutral layer is independent of origin. For a 
supercell containing charged layers like those associated 
with an I-V material, this sum needs to be modified to 
{^jP^j + crjZQj)/c in order to remain meaningful as a 
total polarization. Here (Jj — ±e/S for a Nb02 or KO 
layer respectively, or for a Ti02 or SrO layer. How- 
ever, the second term does not take the form of a sum 
over layer contributions. To cast the sum into this form, 
we define a layer-offcentering polarization that reflects 
the displacement of the layer charge aj from the average 
position of its neighbors according to 

Pj' ^ '^J 2 ~ 4 ^^'■''"^ ~ 4 ^OJ+i ) ■ 

In an extended region of I-V layers, the sum J^j p]° 
counts each charged layer once and only once, but this 
counting is violated at the interfaces with the II-IV lay- 
ers. For example, at a Nb02/SrO interface, we have not 
accounted for a charge density of 6/45 in the last Nb02 
layer, and we have incorrectly assigned a charge of —e/AS 
to the first SrO layer. The missing charge is equivalent to 
a charge density of e/2S located midway between the zqj 
values of these two layers. Similarly, a Ti02/K0 inter- 
face is assigned a charge —e/2S located halfway between 
these layers. Thus, the total polarization of the supercell 
can finally be written as + 4"*<')/c, 

where z™^ and crjf* are the positions and charges of the 
extra interface charges. 

The charges ct™* are, of course, nothing other than 
the compositional polar discontinuity charges CTcomp dis- 
cussed earlier. Having accounted for these, we are left 
with total layer polarizations pj = p^j + associated 
with each layer, which thus provide a layer- by-layer pic- 
ture of the polarization in this type of system. To con- 
vert these into local polarizations Pj having units of 
charge per unit area, we let Pj = Pj/cj where Cj = 
{zo,j+i — zo,j-i)/2. 

The results for n = 4 are shown in Fig. 21 These are 
very indicative of the results obtained for all the supercell 
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FIG. 4: (Color online) Layer-by-layer local polarization profile 
of a supercell with n = 4 (four SrTiOs units and four KNbOa 
units), with corresponding results for n — 2 also indicated 
near each of the interfaces. 

sizes examined. The results for n — 2 centered on each 
interface are also shown on the corresponding segments 
in the figure, allowing the similarity in the behavior of 
the layer polarization to be clearly seen. A saw-tooth- 
like variation about an approximately constant value is 
evident in examining layers belonging to a single ABO3 
constituent. The polarization of the AO layer is larger in 
magnitude than the BO2 in both the I-V and H-IV ma- 
terials in the system, as has been observed elsewherCfi^ 
However, there is a noticeable modification of the Ti02 
layer that is adjacent to KO at the Ti02/K0 heteroint- 
erface; the polarization of this layer is enhanced, so that 
the saw-tooth behavior is broken. 

We find the average polarization deep in the KNbOa 
region to be about 0.45 C/m^, while in the SrTiOa re- 
gion it is about — 0.08C/m^. Recalling that Pg of bulk 
KNbOs at this in-plane lattice constant was found to be 



0.42,C/m^, we find that our earlier expectations for the 
behavior of the polarization profile, that PKNbOa should 
be enhanced relative to its bulk Pg while PsrTiOa should 
be small and of opposite sign, are clearly confirmed. 

In summary, our calculations on perovskite superlat- 
tices composed of II- IV and I-V components have shown 
that sharp local polarization discontinuities can be sta- 
bilized at the charged heterointerfaces that are intrinsic 
to such superlattices. Focusing on the SrTiOa/KNbOs 
system, our calculations show that while the charge com- 
pensation is not perfect, the superlattice period can be 
increased to rather large dimensions before the remaining 
uncompensated bound charges drive the system metal- 
lic. The same principles should apply to II-IV/IITIII 
perovskite superlattices. By suitable choice of materials, 
and/or by tuning the polarization via the application of 
epitaxial strain, which can have a strong effect on the 
spontaneous polarization,^'^ it should be possible to ob- 
tain even better compensation. The resulting superlat- 
tices have strongly broken inversion symmetry, and may 
be of interest for piezoelectric, pyroelectric, non-linear 
optical, and other applications. Finally, we have also 
demonstrated how a layer-by-layer Wannier analysis may 
be applied to a perovskite system in which non-neutral 
layers are present, as an important step towards a more 
complete understanding of the local behavior of multi- 
layered superlattices. 
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